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Abstract: In relativistic systems at RHIC energies, new deconfinement signatures emerge 
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CO . and indicate an increasingly clear separation between soft hadronic processes, and hard 
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{Sj ■ partonic interactions in a locally deconfined subsystem. Here the emphasis is on longitu- 
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dinal variables, namely, net-baryon rapidity distributions. As described in a Relativistic 
Diffusion Model, they change from bell-shaped at the lower to double-humped at the 
higher SPS-energy, but do not reach local statistical equilibrium. At ^snn = 200 GeV 
in the Au-Au system, however, they are shown to consist of three components. In addi- 
QJ , tion to the nonequilibrium contributions, a third fraction close to midrapidity containing 
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Z eq ~ 22 protons reaches local statistical equilibrium in a discontinuous transition. It 
may be associated with a deconfinement of the participant partons and thus, serve as a 
signature for Quark-Gluon Matter formation. 
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1 Introduction 



Considerable progress has recently been made in the long-standing attempts 1J to recre- 
ate a quark-gluon plasma under laboratory conditions by means of relativistic heavy-ion 
collisions. Quark Matter is believed to constitute an important intermediate state of the 
very early universe. There it is in thermal equilibrium, but expands, cools and ends in the 
most dramatic event of a quark-hadron phase transition - a thermal confinement transi- 
tion - at about 10/is: Quarks and gluons condense to form a gas of nucleons and light 
mesons, the latter decaying afterwards. 

As the temperature drops to about 170 MeV, the hadron gas becomes sufficiently 
dilute, and the hadron abundances for different particle species remain fixed ("chemical 
freezeout"). The phase transition has thus set the stage for the subsequent primordial 
synthesis of light nuclei (d, He-3, He-4, Li- 7) at times t~ls, and temperatures T~l MeV 
in the evolution of the early universe [3] . 

The attempt to investigate the quark-hadron phase transition in reverse order, as a 
deconfinement transition, together with the subsequent confinement transition in the labo- 
ratory by means of relativistic heavy-ion collisions experiences severe difficulties. The sys- 
tem has to be sufficiently extended, small collision partners such as protons or deuterons 
are not suitable. The energy density has to reach values above the critical one (about 1.5 
GeV/fm 3 ) for the transition to occur. 

The typical time scale for relativistic heavy-ion collisions in the laboratory is only 
about 10~ 23 s - to be compared with the much larger time scale of the cosmological 
QCD transition of 10 _5 s. Hence, it can not be expected that thermal equilibrium, 
which governs the physical description of the early universe, remains a valid concept 
for theoretical models of relativistic heavy-ion collisions. Whereas particle abundances 
have been shown to be described rather accurately by phase-space ("thermal") models 
jH El E] this does not necessarily imply that the system has reached, or gone through, 
thermal equilibrium. Instead, one has to look for stages of local kinetic equilibrium in 
the short time evolution of the system, and for the possibility that the deconfinement 
transition occurs in such a stage of local thermal equilibrium, affecting only a relatively 
small number of nucleons in a relatively big system. 
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In the fixed-target experiments at the SPS with heavy systems - in particular, with 
the Pb-Pb system at ^s^n = 17.3 GeV - a number of possible phase-transition signatures 
such as strangeness enhancement, and excess of dileptons with invariant mass below that 
of the p-meson had been discussed. The most promising signal, namely, the suppressed 
production of the J/\l/-meson in the presence of a quark-gluon plasma due to vanishing 
string tension and screening, had been predicted by theorists [2] and identified at the SPS 
in heavy systems, but since it could also be caused by hadronic final state interactions 
(nuclear absorption) it seemed not fully convincing. 

Whether the "extra suppression" which was then detected in the Pb-Pb system at 
y / J/v r /v = 17.3 GeV and which could not easily be accounted for by absorption constitutes a 
qgp-signature ist still a matter of debate. At the RHIC energy of 200 GeV per particle, the 
PHENIX collaboration has presented preliminary results for the J/\l/-meson [2j showing 
a slight suppression. In view of the large error bars, however, this is not yet conclusive 
either, one has to wait for more precise data. 

In a probably more promising effort, the four RHIC Au-Au experiments have carefully 
investigated the particle production in central collisions at high transverse momenta. 
When compared to p-p data that are scaled with the number of binary collisions, a 
significant suppression of the produced hadrons is found, which is interpreted as a final 
state effect of the produced dense medium - and possibly, of a quark-gluon plasma. The 
effect may be due to "jet quenching" : Energetic partons traversing the dense medium lose 
energy or are completely absorbed, and the remaining observed hadronic jets are mostly 
created from partons produced near the surface and directed outwards. 

The effect is not observed (instead, the inclusive yield is slightly enhanced) in the 
lighter d-Au system, where compression and heating is much less pronounced, and where 
qgp-formation is therefore unlikely [HI E3 EJ. Moreover, back-to-back pairs are also 
strongly suppressed in central Au-Au for similar reasons, whereas near-side pairs exhibit 
jet-like correlations that are similar to p-p- and d-Au-results, also pointing towards jet 
absorption in the dense (qgp?) medium [T2*] . 

Whereas these results focus essentially on transverse variables, longitudinal variables 
also offer very interesting conclusions regarding hadronic vs. partonic interactions, and 
the possibility of qgp formation. In this contribution, I shall therefore concentrate on 
longitudinal variables - in particular, rapidity distributions of net baryons. The analysis 



3 



is based on a Relativistic Diffusion Model (RDM) that allows to deal with analytical 
solutions, rather than numerical codes that often provide little insight into the physical 
assumptions. 

2 Relativistic diffusion model 

The Relativistic Diffusion Model emphasizes the nonequilibrium-statistical features of 
relativistic heavy-ion collisions. It also encompasses kinetic (thermal) equilibrium of the 
system for times that are sufficiently larger than the relaxation times of the relevant 
variables such as transverse energy Ej_ or rapidity y. 

A first (linear) version of the RDM had been proposed in 1996 and applied successfully 
to the analysis of AGS- and SPS-data, with an emphasis on transverse energy distributions 
integrated over all particle species, and the mean value of the rapidity as function of 
transverse energy Although transverse energy spectra of produced particles turn 

out to be close to thermal equilibrium, some deviations from equilibrium appear in the 
transverse variables. Based on analytical solutions of a transport equation, accurate 
predictions of transverse energy spectra were made, such as in case of the Pb-Pb system 
at the highest SPS energy. 

Distributions of longitudinal variables are of greater interest in the RDM-approach 
since they remain farther away from thermal equilibrium. This is particularly true for 
the rapidity y=l/2 ln((E+p)/(E-p)), which is the Lorentz-invariant counterpart of the 
velocity in the beam direction at relativistic energies. Hence, I have focussed the RDM 
in 1999 on rapidity distributions [18J. For net (participant) baryons, 5-function initial 
conditions corresponding to the beam rapidities are appropriate, and analytical solutions 
of the rapidity transport equation can be compared with data for net-proton rapidity 
spectra. 

Although the initial conditions are less staightforward for produced hadrons, Biyajima 
et al. have started in 2002 |2Z| to use the analytical RDM - which they had developed 
independently from, but with exactly the same result as in - for produced hadrons. 

Comparing to a vast amount of RHIC data for produced charged hadrons (PHOBOS and 
BRAHMS) at both 130 GeV and 200 GeV center-of-mass energy per particle pair, they 
obtain high-precision fits of the data with adjusted values of the friction coefficient (ra- 
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pidity relaxation time in my terminology) and the variance (rapidity diffusion coefficient) 
[2Z1I2H1- With the RDM- approach, they are also able to relate pseudorapidity and rapidity 
distributions to each other, and to establish scaling of the charged-particle rapidity dis- 
tributions with the number of charged hadrons (rather than the number of participants, 
or binary collisions). 

Based on rapidity diffusion coefficients that are not fitted, but instead calculated 
analytically using a dissipation-fluctuation theorem in the weak-coupling limit (where 
the time between two subsequent interactions is large compared to the duration of an 
individual interaction), I have in 1999 obtained good RDM-results for net proton distri- 
butions at the low SIS energies of about 1 GeV per particle. However, this limit is not 
attained at AGS and SPS energies, where progressively larger deviations between RDM 
weak-coupling result and data occur [19J. These deviations have been confirmed in an 
independent numerical calculation by Lavagno in 2002 21 . 

Two solutions of this problem have been offered so far. In a strong-coupling treat- 
ment, the diffusion coefficient D y in rapidity space becomes time-dependent PHIE]- For 
certain parametrizations of D y (t), analytical solutions of the RDM are still possible. As a 
simple substitute in comparisons with data, the enhancement factor due to multiparticle 
(collective) effects may be determined from the deviation between weak-coupling solution 
and data. Typical results for net-proton rapidity distributions at the lower SPS beam 
momentum of 40 GeV/c, the higher SPS beam momentum of 158 GeV/c per particle, 
and RHIC (100 GeV/c per particle in each beam) are shown in Fig.l. 

Alternatively, one may resort to non-extensive statistics, with the underlying relativis- 
ts diffusion equation in rapidity space becoming nonlinear. This approach has been used 
by Alberico et al. in 2000 for transverse mass spectra and transverse momentum fluctu- 
ations at SPS energies, assuming that they were in statistical equilibrium [321 ■ In view 
of the discrepancy between the nonequilibrium weak-coupling result and the SPS Pb-Pb 
data, the approach has been extended to the nonequilibrium situation in rapidity space, 
where Tsallis' nonextensivity parameter [2E] q has been determined from the 158A GeV 
Pb-Pb data [2*T] . Whereas this may well be a reasonable phenomenological parametriza- 
tion of the data, it appears to cover up the multiparticle effects which emerge explicitely 
in the linear approach when comparing data and weak-coupling dissipation-fluctuation 
theorem. 
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Subsequently both the linear (extensive) and the nonlinear (non-extensive) approach 
to nonequilibrium processes in relativistic many-body systems are pursued further, with 
an emphasis on the most recent high-energy results which are investigated experimentally 
at RHIC. I concentrate here on net baryons, and refer the reader to Biyajima et al. for 
produced hadrons in the linear RDM [23 [2E1 E2| • 

3 Net baryon rapidity spectra 

Rapidity distributions of participant (net) baryons are very sensitive to the dynamical 
and statistical properties of nucleus-nucleus collisions at high energies. Recent results 
for net-proton rapidity spectra in central Au+Au collisions at the highest RHIC energy 
of y / s/v r /v = 200 GeV show an unexpectedly large rapidity density at midrapidity. The 
BRAHMS collaboration finds [14 dN/dy = 7.1 ± 0.7 (stat.) ±1.1 (sys.) at y = 0. 

The A, A feed-down corrections reduce this yield by 17.5 per cent when performed 
in accordance with the PHENIX A— results [IS] at 130 GeV, but the amount of stopping 
remains significant, although a factor of about 4 smaller as compared to Pb-Pb at the 
highest SPS energy. (A corresponding STAR result [IB] for y=0 at 130 GeV does not yet 
include the feed-down correction). Many of the available numerical microscopic models 
encounter difficulties to predict the net-proton yield in the central midrapidity valley of 
the distribution, together with the broad peaks at the detected positions. 

Here I interpret the data in the nonequilibrium-statistical Relativistic Diffusion Model. 
The net baryon rapidity distribution at RHIC energies emerges from a superposition of 
the beam-like nonequilibrium components that are broadened in rapidity space through 
diffusion due to soft (hadronic, low p±) collisions and particle creations, and a statistical 
equilibrium (thermal) component at midrapidity that arises from hard (partonic, high 
P±) processes |3*U] . 

At RHIC energies, the underlying distribution functions turn out to be fairly well 
separated in rapidity space. Since the transverse degrees of freedom are in (or very close 
to) thermal equilibrium, they are expected to decouple from the longitudinal ones. The 
time evolution of the distribution functions is then governed by a Fokker-Planck [T7j 
equation (FPE) in rapidity space [HI EH 123 122 121 123 EH • In the more general case of 
nonextensive (non-additive) statistics [29J that accounts for long-range interactions and 
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violations of Boltzmann's Stofizahlansatz [2UI221E2] as wen as f° r non-Markovian memory 
(strong coupling) effects [3TJ E2] , the resulting FPE for the rapidity distribution function 
R(y,t) in the center-of-mass frame is jH0j 



Of 



[R(y,t)} 




J(y)[R(y,t)r\+D(t) — [R(y,t)r. 



(1) 



Here, the rapidity diffusion coefficient D(t) may in general be time-dependent, although 
I will use a constant D y in most of the applications in this work. It accounts for the 
broadening of the rapidity distributions due to interactions and particle creations, and 
it is related to the drift term J(y) by means of a dissipation-fluctuation theorem which 
will be used to actually calculate D r The drift J(y) determines the shift of the mean 
rapidities towards the central value, and I shall discuss linear and nonlinear forms of this 
drift function. 

In derivations of generalized FPE's from the Boltzmann equation, a nonlinear equation 
{n,v 7^1; nonlinear drift function) could in principle be traced back to the nonlinearities 
in the transition probabilities between single-particle states |2S] • However, this would not 
yet include non-Markovian effects. Instead, (JTJ) offers the possibility to describe strong- 
coupling systems that are beyond the realm of the Boltzmann equation. 

Since the norm of the rapidity distribution has to be conserved, \x = 1 is required 
here. It is convenient to introduce a nonextensivity parameter that governs the shape 
of the power- law equilibrium distribution, q = 2 - v I n statistical equilibrium, 

transverse mass spectra and transverse momentum fluctuations in relativistic systems 
at SPS-energies v /saT/v=17.3 GeV require values of q very slightly above one, typically 
q = 1.038 for produced pions in Pb-Pb [32] • For q — > 1, the equilibrium distribution 
converges to the exponential Boltzmann form, whereas for larger values of q (with q < 
1.5) significantly broader equilibrium distributions are obtained, and the time evolution 
towards them becomes super diffusive 

To study rapidity distributions in multiparticle systems at RHIC energies in a 
nonequilibrium-statistical framework [T%1 IT^l |2*T1 12*21 |2U EHII, I start with q = v — 1 
corresponding to the standard FPE. For a linear drift function 



J(y) = (yeq-y)/r t 



(2) 
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with the rapidity relaxation time r y , this is the so-called Uhlenbeck-Ornstein process 
applied to the relativistic invariant rapidity. The equilibrium value is y eq = in the 
center-of-mass for symmetric systems, whereas y eq is calculated from the given masses 
and momenta for asymmetric systems. Using 5— function initial conditions at the beam 
rapidities ±yj, (±5.36 at p=100 GeV/c per nucleon), the equation has analytical Gaussian 
solutions. The mean values shift in time towards the equilibrium value according to 

< Vi,2(t) >= Veq[l ~ exp(-2t/r y )] ± y b exp (-t/r y ). (3) 

For a constant diffusion coefficient D y , the variances of both distributions have the well- 
known simple form 

°Ut) = D y r y [l-ex V (-2t/T y )], (4) 

whereas for a time dependent diffusion coefficient D y {t) that accounts for collective (mul- 
tiparticle) and memory effects the analytical expression for the variances becomes more 
involved [3*T] . At short times t/r y « 1, a statistical description is of limited validity due 
to the small number of interactions. A kinematical cutoff prevents the diffusion into the 
unphysical region \y\ > y b . For larger values of t/r y , the system comes closer to statistical 
equilibrium such that the FPE is valid. 

Since the equation is linear, a superposition of the distribution functions JH] emerging 
from Ri, 2 {y, t = 0) = 5(y T Vb) 



{y-<y x {t)>)\ , 1 r (y-<y 2 (t)>) 2 - 



-exp 



2a 



2 



R(y, t) = — F exp 

2/(2naj(t)) 1 J 2^2^)) 

(5) 

yields the exact solution (normalized to 1). With the the total number of net baryons (or 
protons, depending on the experimental results) Ni+N 2 the rapidity density distribution 
of net baryons (protons) becomes 

dN(y,t = r mt ) = N + N 

dy 

The value of t/r y at the interaction time t=Ti nt (the final time in the integration of (jTj)) 
is determined by the peak positions JH] of the experimental distributions. The same 
approach has also been applied successfully by Biyajima et al. to produced particles at 
RHIC energies |27 3 , although the initial conditions are less straightforward for produced 
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particles, as compared to participant baryons (they use 5-function initial conditions also 
for produced hadrons). 

The microscopic physics is contained in the diffusion coefficient. Macroscopically, the 
transport coefficients are related to each other in the weak-coupling limit (D^ 1 ) through 
the dissipation-fluctuation theorem (Einstein relation) with the equilibrium temperature 
T 

D™ = a-T~f(r y ,T). (7) 

In [TH] I have obtained the analytical result for D™ as function of T and r y from the 
condition that the stationary solution of (0) is equated with a Gaussian approximation to 
the thermal equilibrium distribution in y-space (which is not exactly Gaussian, but very 
close to it) as 



I 



T _T 



D >^ T ) = w C W> T ) ■ (! + 2 - + 2 (-) 2 ) ' ™P( 



m m 



-2 ,2m. 



T 



2^ 

with C(y/s,T) given in m closed form: For a given equilibrium distribution, the 
rapidity diffusion coefficient is then determined by the dissipative constant r y . The proton 
mass is m, and C(y/s,T) ensures that the corresponding thermal distribution function 
Rth is normalized to 1 for each temperature T 

+O0 

R th {y)dy = 1. (9) 



This yields 

+°o i \ ^ 

C{^T)=y b [J (l + 2 XT ( 2/ ) + 2 XT ( 2/ ) 2 )exp(^— ) dy]'" (10) 



with the beam rapidity y& in the center-of-mass, and 



* T M = ~T7~T ( n ) 

(Note that in ^H] the result is written for the beam rapidity in the laboratory system 
yi=5.83 at the SPS-energy for Pb-Pb, corresponding here to y& = ±2.195 in the center- 
of-mass) . 

At fixed beam rapidity the diffusion coefficient displays the expected behavior, 
namely, it rises with increasing temperature of the corresponding equilibrium distribution 
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as in the Einstein relation; the rise is almost linear with T. This allows to maintain the 
linearity of the model and hence, to solve the FPE analytically, although small corrections 
are to be expected. They cause minor deviations in the calculated rapidity distributions 
that are within the size of the error bars of the experimental data at SPS energies. 



4 Strong-coupling diffusion: Memory effects 

In the linear model, net baryon rapidity spectra at low SIS-energies (about 1 GeV per 
particle) are well reproduced, whereas at AGS and SPS energies I find discrepancies jTH] 
to the data that rise strongly with yfs. The origin are most likely strong-coupling effects 
at high energy: the time between two subsequent interactions becomes smaller than the 
duration of an individual interaction, such that the system becomes non-Markovian and 
develops memory-effects. In a schematic approach that serves to outline the effect analyt- 
ically, one may account for the system memory by a time- dependent diffusion coefficient 
through a relaxation ansatz 

(12) 

Here the diffusion coefficient for weak coupling D™ is approached for times t ^> t s larger 
than the time r s that is characteristic for strong coupling - when all secondary particles 
have been created. D™ is well defined in terms of the temperature of the corresponding 
equilibrium distribution, and the particle mass as discussed above. However, for short 
times in the initial phase of the collision before and during particle production, the strong- 
coupling diffusion coefficient D s y dominates and enhances the diffusion in y-space beyond 
the weak-coupling value. This enhancement rises strongly with incident energy as shown 
in ^Hj- It is decisive for a proper representation of the available data for relativistic 
heavy-ion collisions at and beyond SPS energies. The relaxation equation for the rapidity 
diffusion coefficient with the strong-coupling value as an initial condition can be solved 
as 

D y (t) = D™[1 - exp(--)} + D s y exp(--). (13) 

A lower limit for the characteristic time for strong coupling t s can be estimated from 
the time delay of particle production due to quantum coherence [20 which yields about 
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Dy{t) - Dl 
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0.4 fm/c at SPS energies, or simply from the uncertainty principle, which gives about 
the same value for a pion. The strong-coupling time r s is, however, much larger since 
it refers to many particles that are created from the available relativistic energy. The 
transverse energy relaxation time ^Hl - which is of the order of, but slightly smaller than 
the interaction time at SPS energies - gives a reasonable estimate for the strong-coupling 
time such that 

T s < T int < T y (14) 

at SPS-energies. Not much is presently known about the strong-coupling diffusion coeffi- 
cient Dy except that it is significantly larger than its weak-coupling counterpart at these 
energies. Unlike D y , it can not be derived from the known equilibrium distribution. More- 
over, energy and momentum-conservation are not fulfilled in the strong-coupling phase of 
the collision. Leaving a microscopic model to the future, I presently adjust in the linear 
model T) y (t = r int ) to the data as outlined in the next section. Due to the strong coupling 
at short times t/r y <C 1, D y is initially far above the weak-coupling value D y (J5|L and 
it remains there throughout the interaction time at SPS energies, even though it drops 
cslightly until freezeout occurs. The corresponding variance in rapidity space is obtained 
analytically through a solution of the corresponding differential equation for the variances 
with the time-dependent diffusion coefficient (|13j) contained in the inhomogeneity 

d 



f t a 2 y (t) + 2a 2 y (t) = 2D y (t)r y . 
It can be written as a sum of weak- and strong-coupling contributions, respectively 



(15) 



° 2 y (t) = al(t) + al(t) 



(16) 



with 



= D^r y exp(-2t/T y ){[exp(2t/r y ) - 1 

and 

a 2 s (t) = D*T y exp(-2*/V r 



2r, 



Ty 2T S 



Ty 2t s 



exp(—t/r y 



Ty 2T S . 



exp(—t/r r 



)-l]} (17) 



In the limit of r s — > the fluctuations due to strong coupling vanish, and the remaining 
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weak-coupling result for the variance of the rapidity distribution function attains the 
familiar form (j3J) 

for r s ^0. (19) 

Results for the standard deviations in Y-space (Y=y/y&) of the superposed FPE-solutions 
(JSJ) that build up the baryon rapidity distributions in the relativistic system 158 GeV/c 
per nucleon Pb-Pb are shown in Fig.2. The strong-coupling value (upper curve) would 
persist if no particles were created. Due to particle creation, the actual fluctuations 
pass a maximum and gradually approach the weak-coupling result for large times. The 
interaction time (here for a central collision) corresponds approximately to the maximum: 
freezeout occurs for large fluctuations in rapidity at SPS energies and beyond. 

5 From SPS to RHIC energies 

At SPS energies fSJ, t ne discrepancy between weak-coupling result and data has recently 
been confirmed independently in a numerical calculation j2U E2j based on a nonlinear 
drift 

J(y) — ~ a ' m±sinh(y) = —a ■ p\\ (20) 

with the transverse mass m± = \Jm? + p\, and the longitudinal momentum p\\ . Together 
with the dissipation-fluctuation theorem (0), this yields exactly the Boltzmann distribu- 
tion as the stationary solution of for v — q — 1. The corresponding numerical solution 
with 5— function initial conditions at the beam rapidities is, however, only approximately 
correct since the superposition principle is not strictly valid for a nonlinear drift. Still, the 
numerical result shows almost the same large discrepancy between data and theoretical 
rapidity distribution as the linear model. In a q=l framework, the net proton distribution 
in Pb-Pb at the highest SPS energy requires a rapidity width coefficient J ' D y T y that is 
enhanced beyond the theoretical value (JZJ) by a factor of g(\^s) ~ 2.6 due to memory and 
collective effects [IHJ EEH1 EH , Fig.l. 

Alternatively, a transition to nonextensive statistics [2HJ E21 El] maintaining the weak- 
coupling diffusion coefficient from (|7J) requires a value of q that is significantly larger than 



D w r 



1 — exp(—2t/r y ) 
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one. In an approximate numerical solution of (JTJ with the nonlinear drift (j2*Uj) . q=1.25 has 
been determined for the net-proton rapidity distribution in Pb-Pb collisions at the SPS 
[2UI22]- The only free parameter is q, whereas in the linear q=l case the strong-coupling 
enhancement of D y beyond (JZJ) is the only parameter. 

This value of q in the nonlinear model is considerably larger than the result 
( l{^/sNN) = l■12 extrapolated from Wilk et al. [24 at the SPS-energy ^s^n = 17.3 
GeV. Here, the relativistic diffusion approach is applied to produced particles in proton- 
antiproton collisions in the energy range y/s = 53 GeV - 1800 GeV, and used to predict 
LHC-results. The nonlinearity q > 1 appears to be an essential feature of the pp data. 
The larger value of q in heavy systems as compared to pp at the same NN-center-of-mass 
energy emphasizes the increasing superdiffusive effect of multiparticle collisions both be- 
tween participants, and between participants and produced particles. It is, however, 
conceivable that both a violation of (J7J) due to memory effects, and q > 1 have to be 
considered in a complete description. 

The Au+Au system at RHIC energies is then investigated first in the linear q=l 
model for central collisions (10 per cent of the cross-section). Based on the experience 
at AGS and SPS energies [T8J [THJ EH] > it is expected that the nonequilibrium net-proton 
rapidity spectrum calculated with the weak-coupling dissipation-fluctuation theorem (0) 
underpredicts the widths of the nonequilibrium fractions of the experimental distribution 
significantly. At RHIC energies, the precise value of the enhancement due to multiparticle 
effects remains somewhat uncertain at present since the largest-rapidity experimental 
points are on the edge of the nonequilibrium distribution [T41E2]- Typical solutions of 
the linear FPE for various values of t/r y representing the diffusive time evolution of the 
baryonic system due to interactions and particle creations are shown in Fig. 3. 

In the comparison with the BRAHMS data shown in Fig. 4 (bottom), the temper- 
ature T=170 MeV is taken from a thermal fit of charged antiparticle-to-particle ratios in 
the Au+Au system at 200 GeV per nucleon EE] , and the theoretical value of the rapid- 
ity width coefficient calculated from the analytical expression (6) is ^ D y r y = 7.6 • 10~ 2 . 
The weak-coupling nonequilibrium distributions without enhancement due to multiparti- 
cle effects (dotted curves in Fig. 4) are by far too narrow to represent the data, justifying 
the term "anomalous" for the experimental result. 

The calculated distributions become even slightly narrower when T is lowered in order 
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to account for the fact that the equilibrium temperature in the diffusion model should be 
associated with the kinetic freeze-out temperature (which is not yet precisely known at 
RHIC), rather than the chemical freeze-out temperature. As has been discussed in [IT?] , 
the weak-coupling rapidity diffusion coefficient is proportional to the temperature as in 
the theory of Brownian motion j2H|, D y oc T/r y . Hence, lowering the temperature by 
40 MeV reduces the widths of the nonequilibrium distributions by 12 per cent - which is 
hardly visible on the scale of Fig. 4. 

It was shown in [21] for SPS results that the discrepancy between nonequilibrium 
weak-coupling result and data persists in case of the nonlinear drift (|2Uj) that yields the 
exact Boltzmann-Gibbs equilibrium solution for q=l 

d 3 N d 3 N 

d 3 p dy ■ m^dm^dcf) P ( ) 

Hence, it is expected that the nonlinear drift (|2~Uj) does not improve the situation in the 
q=l case at RHIC energies. 

Instead, an enhancement of the weak-coupling rapidity width coefficient by a factor 
of g(\fs) — 3.7 due to collective and memory effects in the system corresponding to 
a violation of (JJJ) yields a good reproduction of the nonequilibrium contributions with 
Tmt/Vj/=0.26 . However, the midrapidity valley that is present in the data is completely 
absent in the extensive nonequilibrium q=l case, solid curves in Fig. 4 (bottom). 

This remains true in the nonextensive case (1 < q < 1.5), with an approximate 
distribution function [2H1 EH EH 1221 12H HU] that is given by a linear superposition of 
power-law solutions of ((H) 

Ri, 2 (y,t) = [1 - (1 - q)^cosh(y- < y lj2 (t) (22) 

The dashed curves in Fig. 4 (bottom) show the result for q=1.4, T=170 MeV and a mean 
transverse mass < m± >=1.2 GeV. This solution is far from the nonextensive equilibrium 
distribution which would be reached for < 2/ 12 — » oo) >= y eq , and it is significantly 
below the midrapidity data. The result is even worse for larger values of m±. In contrast, 
the Pb-Pb data at SPS energies [10] are well described both in the linear model ^H] (Fig.6, 
top) and in the nonlinear case (cf. [21] for results with a time- dependent temperature 
and an integration over transverse mass). 
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6 Local deconfinement at RHIC 



It turns out, however, that the RHIC data can be interpreted rather precisely in the 
linear q=l framework with the conjecture that a fraction of Z eq ~ 22 net protons near 
midrapidity reaches local statistical equilibrium in the longitudinal degrees of freedom, 
Figs.4(top),5. The variance of the equilibrium distribution R eq (y) at midrapidity is broad- 
ened as compared to the Boltzmann result due to collective (multiparticle) effects by the 
same factor that enhances the theoretical weak-coupling diffusion coefficient derived from 
(J7J), shaded areas in Figs. 4, 5. This may correspond to a longitudinal expansion (flow) 
velocity of the locally equilibrated subsystem as accounted for in hydrodynamical de- 
scriptions. In the nonextensive model, the corresponding local equilibrium distribution is 
broadened (blue-shifted) according to q ~ 1.4. Here the enhanced value of q appears as 
a convenient parametrization of collective expansion ("longitudinal flow"). 

Microscopically, the baryon transport over 4-5 units of rapidity to the equilibrated 
midrapidity region is not only due to hard processes acting on single valence (di) quarks 
that are described by perturbative QCD, since this yields insufficient stopping. Instead, 
additional processes such as the nonperturbative gluon junction mechanism jSJj are nec- 
essary to produce the observed central valley. This may lead to substantial stopping 
even at LHC energies where the separation of nonequilibrium and equilibrium net baryon 
fractions in rapidity space is expected to be even better than at RHIC. In the late ther- 
malization stage [HE], nonperturbative approaches to QCD thermodynamics are expected 
to be important. 

Recent work indicates that one may account for the observed stopping in heavy- 
ion collisions at SPS and RHIC energies with string-model parameters determined from 
hadron-hadron collisions If this was confirmed, the corresponding rapidity distribu- 
tions would not be considered to be anomalous from a microscopic point of view. However, 
this view does not offer a distinction between nonequilibrium and equilibrium contribu- 
tions to the net baryon rapidity spectra, which both exist at RHIC energies, and are 
anomalously broadened. 

Macroscopically, the complete solution of (JTJ) in the q=l case is a linear superposition 
of nonequilibrium and local equilibrium distributions (Fig. 6, bottom) 

R{y, t = T int ) = R^y, r int ) + R 2 {y, r int ) + R l ° c (y). (23) 
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with the same enhancement factor g(y/s) due to multiparticle (collective) effects for all 
three distributions. The net-baryon rapidity distribution becomes 

dN{V,t d y Tmt) = NiRi(y, T int ) + N 2 R 2 (y, r mt ) + N eq R% c (y) (24) 

with the number N eq of net baryons (here: net protons) in local equilibrium near mid- 
rapidity, and Ni + N 2 + N e(? equal to the total number of net baryons (158 net protons for 
central Au-Au). This yields a good representation of the preliminary BRAHMS data [14 . 
(In the q> 1 case, the corresponding solution is questionable because the superposition 
principle is violated). Based on (j24jl . the transition from net-proton rapidity spectra with 
a central plateau in Pb-Pb at the lower SPS energies [3T], via a double-humped distribu- 
tion at the maximum SPS energy gU EH |2T1 122 E3 to the central valley at RHIC [THEOj 
is well understood. It has not yet been possible to identify a locally equilibrated subsys- 
tem of net baryons at midrapidity below RHIC-energies, although it cannot be excluded 
that it exists. At SPS energies, the data jlU] are well described by the nonequilibrium 
distributions, and it is much more difficult (and probably impossible) to identify a locally 
equilibrated component because the relevant rapidity region is comparatively small, and 
an equilibrated contribution cannot be separated from the nonequilibrium components in 
rapidity space. In pp-collisions at yfs =53-900 GeV, no convincing signatures of a phase 
transition were found [4*T] . 

Most remarkably, Figs. 4-6 suggest that in central Au+Au collisions at v /J/v r /v= 200 
GeV there is no continuous transition from the nonequilibrium to the equilibrium contri- 
bution in net-proton rapidity spectra as function of time. The central valley in net-proton 
rapidity spectra at RHIC energies could thus be used as an indicator for partonic processes 
that lead to a baryon transfer over more than 4 units of rapidity, and for quark-gluon 
plasma formation. The discontinuity may well be due to a sudden enhancement in the 
number of degrees of freedom as encountered in the deconfinement of participant partons, 
which enforces a very rapid local equilibration in a fraction of the system as indicated in 
Figs. 4, 5 through a discontinuity in the time evolution of the solutions of the FPE, which 
are afterwards very close to the equilibrium result. Here, the sudden enhancement in 
the number of degrees of freedom by a factor of about 6 in the quark-gluon plasma as 
compared to the hadronic phase [3] is modelled by a corresponding increase in time, Fig.5. 

The linear RDM with three sources j30j including a local equilibrium fraction near 
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midrapidity has recently also been adopted by Biyajima et al. [12] for produced particles 
in analyses of BRAHMS jS] and PHOBOS g3| pseudorapidity data at y/s =130 GeV and 
200 GeV, using the proper Jacobian to transform from y to r]=-lntan(0/2). In particular, 
the PHOBOS data at 130 (200) GeV and 0-6 per cent centrality are well reproduced 
with 3134 (3858) charged hadrons in the central source indicative of an equilibrated qgp, 
whereas only 896 (1102) charged hadrons reside in the non-equilibrated fragmentation 
regions. 

These very recent results for produced charged hadrons are, however, still somewhat 
ambiguous in view of (1) the uncertainties regarding the initial conditions for produced 
particles in rapidity space, and (2) the difficulties to clearly separate fragmentation and 
central regions for produced particles. Nevertheless, the RDM-analyses of produced 
hadrons by Biyajima et al. ^2] provide additional evidence for local equilibrium near 
midrapidity. 

7 Conclusion 

To conclude, I have interpreted recent results for central Au-Au collisions at RHIC ener- 
gies in a Relativistic Diffusion Model (RDM) for multiparticle interactions based on the 
interplay of nonequilibrium and local equilibrium ("thermal") solutions. In the linear ver- 
sion of the model, analytical results for the rapidity distribution of net protons in central 
collisions have been obtained. The anomalous enhancement of the diffusion in rapid- 
ity space as compared to the expectation from the weak-coupling dissipation-fluctuation 
theorem due to high-energy multiparticle effects has been discussed using extensive and 
nonextensive statistics. 

A significant fraction of about 14 per cent of the net protons reaches local statistical 
equilibrium in a fast and discontinuous transition which is likely to indicate parton decon- 
finement. The precise amount of protons in local equilibrium is related to the experimental 
value of the rapidity density close to y=0 and hence, possible changes in the final data 
will affect the percentage. It has not yet been possible to isolate a corresponding fraction 
of longitudinally equilibrated net protons in the Pb-Pb system at SPS energies. Since no 
signatures of a transition to the quark-gluon plasma have been observed in pp-collisions, 
Quark-Matter formation is clearly a genuine many-body effect occuring only in heavy 
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systems at sufficiently high energy density. Consequently, a detailed investigation of the 
flat midrapidity valley found at RHIC, and of its energy dependence is very promising. 
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Figure captions 



FIG. 1: Net proton rapidity spectra in the Relativistic Diffusion Model (RDM), solid 
curves, at the lower SPS momentum of 40 GeV/c per particle (top; data not yet 
available), at 158 GeV/c (compared to NA 49 data [H]|; cf. Fig.6 for error bars) and 
at the highest RHIC energy ^snn = 200 GeV (compared to preliminary BRAHMS 
data Jl]). The variable Y=y/y\ ) renormalizes the y-distributions in the center-of- 
mass with the beam rapidities ±y&. Dashed lines are thermal equilibrium results 
without expansion. The transition from bell-shaped to double-humped is clearly 
shown in the RDM. 

Fig. 2. Standard deviations oy(t) of the solutions of the Fokker-Planck equation that 
build up the Pb-Pb baryon rapidity distribution at SPS momenta of 158 GeV/c 
per particle in Y-space (Fig.l). The upper curve is the strong-coupling result. Due 
to particle creation, the actual fluctuations (middle curve) gradually approach the 
weak-coupling result, lower curve. The interaction time T int /r y (arrow) approxi- 
mately corresponds to the maximum of the fluctuations in rapidity space. 

FIG. 3: Analytical solutions of the Relativistic Diffusion Model for various values of t/r y 
representing the diffusive time evolution of the baryonic system due to interactions 
and particle creations for central Au-Au at RHIC energies. Here the net proton 
content is 136 protons. This subsystem remains in a non-equilibrium state in the 
actual experiment, cf. Fig.6. 

FIG. 4: Nonequilibrium contributions to the net-proton rapidity spectra of Au-Au at 
v /sjv r /v = 200 GeV in the Relativistic Diffusion Model (RDM) with an equilibrium 
temperature of T=170 MeV (bottom). The solid curve is obtained in the linear 
model (q=l). Its width is ("anomalously") enhanced as compared to the theoretical 
weak-coupling value (dotted curves) by g(y/s) = 3.7 due to multiparticle effects 
[T9l I3T] according to the preliminary BRAHMS data [H] , squares. The net proton 
content is 158. The dashed curve corresponds to q=1.4 and < m± >=1.2 GeV in 
the nonlinear model. Reasonable agreement with the midrapidity data requires that 
14 per cent of the net protons reach local thermal equilibrium in a discontinuous 



22 



transition, top (dashed curve: thermal equilibrium result for 22 protons; shaded 
area: broadening due to multiparticle effects with g(^/s) = 3.7). 

FIG. 5: Analytical solutions of the Relativistic Diffusion Model for various values of t/r y 
representing the time evolution of a local subsystem of the net baryons in central 
Au-Au at RHIC energies with a sudden enhancement of the number of degrees of 
freedom due to deconfinement, and the subsequent fast local equilibration. The net 
proton content of this local subsystem is 22. The shaded area is the local equilibrium 
distribution centered at midrapidity. 

FIG. 6: Net-proton rapidity spectra for central collisions of Au-Au at ^s^n — 200 GeV 
consist of two nonequilibrium components (solid peaks, bottom; 136 protons) plus 
a local equilibrium contribution at T=170 MeV. The shaded area shows its broad- 
ening due to collective (multiparticle) effects by the same factor g(\J~s) = 3.7 as 
the nonequilibrium fractions. After hadronization, it contains Z eq ~ 22 protons. 
Superposition creates the flat valley near midrapidity (bottom) in agreement with 
the preliminary BRAHMS data points jHj. Diamonds include A feed-down correc- 
tions at y=0 (17.5 per cent) and y=2.9 (20 per cent), respectively. Arrows indicate 
the beam rapidities ±y b . At SPS energies (top), NA49 data [40J for central Pb-Pb 
events (5 per cent) including A feed-down corrections are compared with the pure 
nonequilibrium result of the linear model fHl E] ■ Here, the net proton content is 
164. 
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